void InitWatAngCoords( FILE *fPtr )
{
  VecR e;
  real eAng[ 3 ];
  //, phi, theta, psi, cphi, ctheta, cpsi, sphi, stheta, spsi;
  int n;

  fprintf( fPtr, "InitWatAngCoords... "); fflush(fPtr);

  DO_SITE {
    if ( site[ n ].mechType == NSRB ) {
      VRand( &e );
      eAng[ 0 ] = atan2( e.x, e.y ); // phi
      eAng[ 1 ] = acos( e.z );       // theta
      eAng[ 2 ] = 2. * PI * RandR(); // psi
      //      EulerToQuat( &site[ n ].q, eAng );
      /*       printf("*** w:% f x:% f y:% f z:% f ok\n", site[n].q.u4, site[n].q.u1, site[n].q.u2, site[n].q.u3 ); */
      EulerToRotMat( &siteRotationMatrix[ n ], eAng[ 0 ], eAng[ 1 ], eAng[ 2 ] ); 
      //      PrintMatrix(site[ n ].qMat);
      /*       Quat qCheck; */
      /*       RotMatToQuat( &qCheck, &siteRotationMatrix[ n ] ); */
      /*       printf("*** w:% f x:% f y:% f z:% f old\n", qCheck.u4, qCheck.u1, qCheck.u2, qCheck.u3 ); */
      /*       RMat checkMat; */
      /*       BuildRotMatrixFromQuat( &checkMat, &qCheck, 1 ); */
      /*       printf("% 4.3f % 4.3f % 4.3f\n% 4.3f % 4.3f % 4.3f\n\n", */
      /* 	  siteRotationMatrix[ n ].u[ 0 ], siteRotationMatrix[ n ].u[ 1 ], siteRotationMatrix[ n ].u[ 2 ], */
      /* 	  checkMat.u[ 0 ], checkMat.u[ 1 ], checkMat.u[ 2 ] ); */
      /*       BuildQuatFromRotMat( &qCheck, &site[ n ].qMat ); // 27/3 trying to fix this... */
      /*       printf("*** w:% f x:% f y:% f z:% f today \n\n", qCheck.u4, qCheck.u1, qCheck.u2, qCheck.u3 ); */
    }
  }
  fprintf( fPtr, "done\n" ); fflush(fPtr);
}

